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We investigate the rates of production and thermalization of and v T neutrinos at temperatures 
and densities relevant to core-collapse supernovae and protoneutron stars. Included are contribu- 
tions from electron scattering, electron-positron annihilation, nucleon-nucleon bremsstrahlung, and 
nucleon scattering. For the scattering processes, in order to incorporate the full scattering kine- 
matics at arbitrary degeneracy, the structure function formalism developed by Reddy et al. (1998) 
and Burrows and Sawyer (1998) is employed. Furthermore, we derive formulae for the total and 
differential rates of nucleon-nucleon bremsstrahlung for arbitrary nucleon degeneracy in asymmetric 
matter. We find that electron scattering dominates nucleon scattering as a thermalization process 
at low neutrino energies (e„ < 10 MeV), but that nucleon scattering is always faster than or compa- 
rable to electron scattering above e„ ~ 10 MeV. In addition, for p > 10 13 g cm -3 , T < 14 MeV, and 
neutrino energies < 60 MeV, nucleon-nucleon bremsstrahlung always dominates electron-positron 
annihilation as a production mechanism for and u T neutrinos. 
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I. INTRODUCTION 



The cores of protoneutron stars and core-collapse supernovae are characterized by mass densities of order ~ 10 10 — 
10 14 g cm -3 and temperatures that range from ~ 1 to 50 MeV. The matter is composed predominantly of nucleons, 
electrons, positrons, and neutrinos of all species. For v„ and v T types (collectively 'i^s'), which carry away 50—60% 
of the ~ 2 — 3 x 10 53 ergs liberated during collapse and explosion, the prevailing opacity and production processes 
are ^-electron scattering, ^-nucleon scattering, electron-positron annihilation (e + e _ <-> v^v^), and nucleon-nucleon 
bremsstrahlung. While all of these processes contribute for the electron types (z/ e s and u e s), the charged-current 
absorption processes v e n <-> pe~ and D e p <-> ne + dominate their opacity so completely that in this paper we address 
only production and thermalization. 

Supernova theorists had long held |jj that i^-nucleon scattering was unimportant as a mechanism for neutrino 
equilibration. While this process was included as a source of opacity f|J|, it served only to redistribute the neutrinos 
in space, not in energy. In contrast, ^-electron scattering was thought to dominate neutrino thermalization. In 
addition, the only v^v^ pair production mechanisms employed in full supernova calculations were e + e _ <-> v^v^ and 
plasmon decay (~f p i <-*■ v^v^) nucleon-nucleon bremsstrahlung was neglected as a source. Recent developments, 
however, call both these practices into question and motivate a re-evaluation of these opacities in the supernova 
context. Analytic formulae have recently been derived [^-^) which include the full nucleon kinematics and Pauli 
blocking in the final state at arbitrary nucleon degeneracy. These efforts reveal that the average rate of energy 
transfer in i/^-nucleon scattering may surpass previous estimates by an order of magnitude Hence, this 

process may compete with i^-electron scattering as an equilibration mechanism. Similarly, estimates for the total 
nucleon-nucleon bremsstrahlung rate have been obtained [Q,|l2]-[l4| which indicate that this process might compete 
with e + e~ annihilation. 
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These results suggest that the time is ripe for a technical study of the relative importance of each process for 
production or thermalization. To conduct such a study, we consider neutrinos in an isotropic homogeneous 
thermal bath of scatterers and absorbers. In this system, the full transport problem is reduced to an evolution of the 
neutrino distribution function (TJ) in energy space alone. Although this is a simplification of the true problem, it 
provides a theoretical laboratory in which to analyze the rates both for equilibration of an initial neutrino distribution 
function with dense nuclear matter and for production of the neutrinos themselves. From these rates we determine 
the importance and particular character of each process, and discover in which energy, temperature, or density regime 
each dominates. We employ a general prescription for solving the Boltzmann equation in this system with the 
full energy redistribution collision term. We compare quantitatively, via direct numerical evolution of an arbitrary 
neutrino distribution function, the rates for thermalization and production by each process, at all neutrino energies. 
Furthermore, we present the total nucleon-nucleon bremsstrahlung rate for arbitrary nucleon degeneracy and derive 
the single and production spectra. This facilitates a more comprehensive evaluation of its relative importance 
in neutrino production than has previously been possible. 

In §||, we discuss t he g eneral form of the Boltzmann equation and our use of it to study equilibration and 



production rates. In § yll , we provide formulae for each of the four processes we consider: z-^-nucleon scattering, 
^-electron scattering, and v^v^ pair production via both nucleon-nucleon bremsstrahlung and e + e~ annihilation. In 
§IV , we present the results of our equilibration calculations, showing the time evolution of distribution functions 
as influenced by each of these processes individually. We include plots of thermalization and production rates for 
each process as a function of neutrino energy and time. For the scattering interactions we include figures of the time 
evolution of the net energy transfer to the medium as a function of incident neutrino energy. We repeat this analysis 
at points in temperature, density, and composition space relevant to supernovae and protoneutron stars, taken from 
snapshots of a stellar profile during a realistic collapse calculation 0. Using these results, we discuss the relative 
importance of each process in shaping the emergent spectrum. In §M, we recapitulate our findings and conclusions. 



II. THE BOLTZMANN EQUATION 

The static (velocity = 0) Boltzmann equation for the evolution of the neutrino distribution function {T v ) 1 including 
Pauli blocking in the final state, and for a spherical geometry, is 

\m + 4r + ^1) = (1 - ™» - (2 - 1} 

where t is the time, r is the radial coordinate, and fi(= cos (9) is the cosine of the zenith angle. j v and Xv are the total 
source and sink, respectively. For emission and absorption, j v is the emissivity and \ v is the extinction coefficient. 
For scattering, both j„ and Xu are energy redistribution integrals which couple one neutrino energy bin with all 
the others. The matrix element and associated phase-space integrations which comprise j v and Xv for electron and 
nucleon scattering yield the probability that a given collision will scatter a particle into any angle or energy bin. A 



full transport calculation couples energy and angular bins to each other through the right hand side of eq. (2.1). 

In a homogeneous, isotropic thermal bath of scatterers and absorbers no spatial or angular gradients exist. Conse- 
quently, the Boltzmann equation becomes 

-°4^ = {\-T v)3v -T v x v - (2.2) 
c at 

By dealing with this system, the transport problem reduces to an evolution of T v in just energy and time. Note 
that for scattering processes, both j v and %„ require an integral over the scattered neutrino distribution function 
T' v . Similarly, in evolving T v via the production and absorption processes, j v and Xv involve an integration over the 
anti- neutrino distribution function !F 9 . Therefore, T„ must be evolved simultaneously with T v . While j v and Xv may 



be fairly complicated integrals over phase-space, the numerical solution of eq. (2.2) is straightforward. 

Given an arbitrary initial !F V , we divide the relevant energy range into n energy bins. We then solve eq. (2.2) for 
each bin individually and explicitly. Angular integrals over scattering cosines which appear in the f^-nucleon and 
^-electron scattering formalism, as well as the electron energy integration needed for e + e~ annihilation, are carried 
out with a 4-point Gauss-Legcndre integration scheme. The double integral over dimensionless nucleon momentum 
variables needed to obtain the contribution from nucleon-nucleon bremsstrahlung is computed using nested 16-point 
Gauss-Laguerre quadratures. 
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A. Rates for T v Evolution and Energy Transfer 



Scattering, emission, and absorption processes, at a given neutrino energy (e„), produce and remove neutrinos 
from the phase-space density at that energy. The former achieves this by transferring energy to the matter during 
scattering, the latter two by emitting or absorbing directly from that bin. The Boltzmann equation can then be 
written in terms of an in and an out channel, the former a source and the latter a sink 



dt 



dT v 



Of 



dT v 



dt 



(2.3) 



Consequently, for any interaction, there are two rates to consider: the rate for scattering or production into a given 
energy bin (Ti n ) and the inverse rate for scattering or absorption out of that bin (T ou t). The rates cj v and cx„ yield 
timescales for an interaction to occur, but fail, in the case of the former, to fold in Pauli blocking in the final state. 
Equation (2.3) includes these effects and provides natural timescales for T v evolution: 



-L in. ^~ 



T v dt 



■ CJu 



and 



1 dT v 



T v dt 



CXu- 



(2.4) 



(2.5) 



Note that although eq. ( |2.5| ) does not explicitly contain a Pauli blocking term, x„ contains an integral over (1 — T'J), 
in the case of scattering, and an appropriate final-state blocking term, in the case of absorption. At a given e v , then, 
Ti n incorporates information about the v u phase-space density at that energy. Conversely, at that same e„, T out 
contains information about the phase-space at all other energies. Regardless of the initial distribution, dT v jdt = in 
equilibrium. This implies I\ n — T out in equilibrium and, hence, we build in a test for the degree to which the system 
has thermalized. 

Just as there are distinct rates for the in and out channels of the Boltzmann equation during equilibration, so too 
are there distinct scattering energy transfers. For scattering with a scatterer s (electron or nucleon), at a specific 
£„, two thermal average energy transfers can be defined; 



and 



(cj) out = / d 3 P ' u Cj(l-F' v )l° Ut [u M 8- 



(2.6) 



(2.7) 



where primes denote the scattered neutrino, u>(= e v — s' u ) is the energy transfer, and X m and T out are the kernels for 
scattering into and out of a given energy bin, respectively. As a consequence of detailed balance between the in and 
out channels of the Boltzmann equation, I™ — e _,3w X otlt , where f3 — l/fc^T and T is the matter temperature. (The 



scattering kernels are discussed in detail in §111 for both scattering processes.) Note that the denominators in eqs. 
(2.6) and (2.7), up to constants which divide out in the definitions of {oS)i n and (u>) ut, are just j v and \v, respectively. 



In an effort to provide more than one measure of the timcscale for T v equilibration due to scattering and to make 
contact with previous neutrino thermalization studies we also define a set of timescales in terms of (ui) ou t and 

the higher w-moment, (uj 2 )out', 



and 



(<jj)out 



(2.8) 



(2.9) 



Td is the rate for shifting the centroid of a given distribution and is the rate for spreading an initial distribution 
p5| . In contrast with the work of Jl5|-[T7|, we include the full effects of Pauli blocking in the final state, allowing us 
to deal consistently with cases in which the v^s are partially degenerate. 
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III. INDIVIDUAL INTERACTIONS 



This section details the source and sink terms necessary to solve the Boltzmann equation for the time-evolution of 
T v . Sections § [III A and §HIB are dedicated to the presentation and discussion of the collision terms for ^-nucleon 
and f^-electron scattering, respectively. Section § [II C describes the Legendre series expansion approximation and the 
use of it to compute the contribution to the Boltzmann equation, the pair emissivity, and the single f M spectrum due to 
e+e~ <-> v^v^. Our derivations of j„ and x„, as well as the single and pair spectra from nucle on-nucleon bremsstrahlung 
at arbitrary nuclcon degeneracy and in the non-degenerate limit, are presented in §111 D. In what follows, we take 



G 2 



1.55 x 10" 



MeV s , sin 6w — 0.231, and employ natural units in which % = c = ks = 1- 



A. Nucleon Scattering: u^n <-> v^n and v^p «-» v^p 

Researchers working on supernova and protoneutron star evolution have recently re-evaluated the issue of energy 
transfer via f M -nucleon scattering |p|,[7|-^o| . Originally, the assumption was made that the nucleons were stationary 
Q. If a neutron of mass m n is at rest with respect to an incoming neutrino of energy s v , one finds that the energy 
transfer (u>) is ~ —e 2 /m n . For e v — 10 MeV, Q ~ —0.1 MeV, a fractional energy lost of 1%. Using these simple 
kinematic arguments and disregarding neutrino and nucleon Pauli blocking, one finds that the thermalization rate 
for f M -electron scattering should be approximately a factor of 20 larger than that for i/^-nucleon scattering. In the 
context of interest, however, at temperatures of order 10 MeV and mass densities of order 10 13 g cm -3 , free nucleons 
are not stationary, but have thermal velocities. The fractional energy exchange per collision, in the case of ^-neutron 
scattering, is then ~ p n /m n c 0. For T ~ 10 MeV this gives a ^10—20% change in e„ per collision. This calls the 
naive estimate of the relative importance of ^-nucleon scattering as a thermalization process into question and a 
more complete exploration of the relative importance of the two scattering processes is necessary. 

Recently, analytic formulae have been derived which include the full kinematics of ^-nucleon scattering at arbitrary 
nucleon degeneracy |4|-|(|. At the temperatures and densities encountered in the supernova context non- interacting 
nucleons are not relativistic. Due to nucleon-nucleon interactions, however, at and around nuclear density (~ 2.68 x 
10 14 g cm -3 ), the nucleon's effective mass drops and is expected to be comparable with its Fermi momentum [EJ. 
In such a circumstance, a relativistic description of the ^ M -nucleon scattering interaction is warranted. In addition, 
spin and density correlation effects engendered by these nucleon-nucleon interactions have been found to suppress the 
j/u-nucleon interaction rate by as much as a factor of ~ 2 — 3 S[LG,[ll| . 



In this study, we focus on equilibration rates at densities < 1 x 10 14 g cm -3 where it is still unclear if nucleon- 
nucleon interactions will play an important role. This ambiguity is due in part to uncertainties both in the nuclear 
equation of state and the nucleon-nucleon interaction itself. For this reason we choose to treat the nucleons as non- 
relativistic and non-interacting, thereby ignoring collective effects which might enhance or reduce the ^-nucleon 



scattering rate. Making these assumptions, we find that j v and Xv m eq. ( |2.2| ) are given by 

G 2 



3* 

and 

Xv - 



(27T) 



G 2 
(2ny 



Jd 3 p^l NC K^ (3-1) 



[ d 3 ftl NC [l-K}, (3-2) 



where f3 — 1/T, p' v is the final state neutrino momentum, and u is the energy transfer. In eqs. (3.1) and (3.2), the 
neutral-current scattering kernel is given by 

X NC = S(q, lo) [(1 + li)V 2 + (3 - n)A 2 ], (3.3) 

where cos9) is the cosine of the scattering angle between incident and final state neutrinos and S(q,u>) is the 



dynamic structure function. In eq. (3.3), V and A are the appropriate vector and axial-vector coupling constants; for 



i/ M -neutron scattering, V = —1/2 and A = —1.26/2. The dynamic structure function is 

S(q, u,) = 2jj^F(l-F') (2^(c +e-e') 



(2tt) 3 

= 2 Iran (0) (q, u) (1 - e-^y 1 (3.4) 
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where q = \p v — p' v 



[el 



,12 



2e,e>] 1 /2 is 

the magnitude of the momentum transfer, and J- and T 1 are the 



incident and scattered nucleon distribution functions, respectively. In eq. (3.4), p is the incident nucleon momentum, 
e is the incident nucleon energy, and e' is the scattered nucleon energy. The imaginary part of the free polarization 
is given by |i 18 



Imn (0) (g,w) 



2n(3q 



-In 



1 



1 -f e -Q 2 +r)-,3u 



(3.5) 



where 



Q 



m(3 



1/2 



U) 

Q 



q 

2m 



(3.6) 



77 is the nucleon degeneracy (/x/T), and m is the nucleon mass. The factor e _/3w which appears in eq. (3.1 
consequence of the fact that S(q, — uo) = e~" u S(q, cu), itself a consequence of detailed balance between the in and 
out channels of the Boltzmann equation. The dynamic structure function can be thought of as a correlation function 
which connects £„ and e' v . 

The <fi angular integrations implicit in eqs. (3.1) and (3.2) can be computed trivially assuming the isotropy of T v . 
Furthermore, defining the coordinate system with the momentum vector of the incident neutrino, the scattering angle 
and the direction cosine are equivalent. Combining these two equations in the Boltzmann equation for the evolution 
of T v due to neutral-current ^-nucleon scattering, we obtain 



G 2 



dt ~ (2tt) 2 



J de' u e' 2 J ^ dfi X NC { [1 - ?v\ K^ (iu> -T v \\- T' v \ } 



(3.7) 



B. Electron Scattering: v^e 



At the temperatures and densities encountered in supernovae and protoneutron stars, electrons are highly rela- 
tivistic. A formalism analogous to that used for ^-nucleon scattering is desired in order to include the full electron 
kinematics at arbitrary electron degenera cy. R eddy et al. Q have developed a relativistic generalization of the 
structure function formalism described in § [II A . They obtain a set of polarization functions which characteriz e th e 
relativistic medium's response to a neutrino probe in terms of polylogarithmic functions. In analogy with eq. (fe.Tf) , 
we can write the Boltzmann equation for the evolution of T v due to i^, -electron scattering, as 



dt 



(3.8) 



where XjJ c is the relativistic neutral-current scattering kernel for v^s, analogous to 2nc in eq. (3.3). All the physics 
of the interaction is contained in T^ c , which can be written as 



I N r c =Im{A^n« }(1 



r 1 . 



As in the non-relativistic case, I^q is composed of the lepton tensor, 

K af3 = 8[2fc Q fc' 3 + (k ■ q)g al3 - (fcV + q a k ) - ie a ^ u k^q v ], 



(3.9) 



(3.10) 



which is just the squared and spin-summed matrix element for the scattering process written in terms of k a , the 
incident four-momentum, and q a (= (u>,q)), the four-momentum transfer. The scattering kernel also contains the 
retarde d po larization tensor, H^g, which is directly analogous to the free polarization in the non-relativistic case given 
in eq. (3.4). The retarded polarization tensor is related to the causal polarization by 



Imll^a = tanh I — ^-(3co I Imll^ 



(3.H) 



and 
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/i4 
—^Tr[G e (p)J a G' e (p + q)J }, (3.12) 

In eq. Q, p a is the electron four-momentum and J a is the current operator. The electron Green's functions (G e 
and G" e ), explicit in the free polarization, connect points in electron energy space and characterize the effect of the 
interaction on relativistic electrons. The polarization tensor can be written in terms of a vector part, an axial-vector 
part, and a mixed part, so that 

n Q/3 = v 2 u^ + A 2 n* fj - 2vau v J. (3.13) 

In turn, the vector part o f th e polarization tensor can be written in terms of two independent components, ELr and 
IT^. In contrast with eq. (3.3), since v/c ~ 1 for the electrons, the angular terms which were dropped from the matrix 
element in the non-relativistic case, leading to a single structure function, must now be retained. can then be 
written as a set of three structure functions Q : 

= 8[AS 1 (q,u;)+S 2 (q,cj) + BS 3 (q,cj)](l - e^")" 1 , (3.14) 

where A — (4e„e[, + q 2 )/2q 2 and B = e v + e' v . These structure functions can be written in terms of the vector parts 
of the retarded polarization tensor (11^ and Il£), the axial part (n^), and the mixed part (Hy A ): 

Si(q,u) = (V 2 +A 2 ) [Imn£(g,w)+Imn£(g,a;)], (3.15) 



S 2 (q,w) = {V 2 +A 2 ) ImITf( 9) u;)- A 2 ImU^{q,uj), (3.16) 

and 

S 3 (q,w) = 2VATmn$ A (q,w). (3.17) 

The retarded polarization functions, in terms of differences between polylogarithmic integrals, can be found in Ap- 
pendix |a|. 



C. Electron- Positron Annihilation: e e <-> v^v^ 

Fermi's Golden Rule for the total volumetric emission rate for the production of z/^s via electron-positron annihi- 
lation can be written as 



where 



E\f\ = (1 - - F P )F e -F e+ , (3.19) 



and S 4 (P) conserves four-momentum. In eq. (3.18), p a {= (e,p)) and p' a {= (e',p')) are the four-momenta of the 
electron and positron, respectively, and <?"(= (e u ,ql)) and q"(= (ep,q„)) are the four-momenta of the and D^, 
respectively. The process of electron-positron annihilation into a neutrino/anti-neutrino pairs is related to neutrino- 



electron scattering considered in §111 B via a crossing symmetry. In order to make the problem tractable, we follow the 
standard procedure Q of expanding the production kernel in a Legendre series in the scattering angle to first order (see 
Appendix ^). Near the neutrinospheres, at densities which render neutrino transport diffusive this approximation 
holds. In a full neutrino transport algorithm, however, which must handle both the diffusion and free-streaming 
limits, the second-order term, with proper closure relations, must be used in the semi-transparent regime between 
the neutrinospheres and the shock Jl9| . Having made this approximation, including only the zeroth- and first-order 
terms, the single spectrum is 

^ = (1 ^ j°° d ^ 4 *o(^) (1 - ft), (3-20) 



G 



where $q(e v ,£„) is the zeroth-order production kernel expansion coefficient, an integral over the electron energy (see 
Appendix |bJ) B. With the differential spectrum or emissivity (dQ/de v ) in hand, it is a simple matter to extract the 
contribution to the Boltzmann equation due to e + e~ annihilation. As eq. ( [3.20 ) already contains the blocking 
factor, the contribution to the Boltzmann equation, the in channel explicit in eq. (|2.3|), can be written as pi 



f (2tt) 3 dQ 
4n e 3 de u ' 



(3.21) 



In order to obtain the out channel for absorption due to e + e~ annihilation, we need only replace T e -J- e + in eq. (3.19) 
with an electron/positron blocking term, (1 — T e -){\ — J~ e +), and replace the and blocking terms in eq. (3.2C) 
with T v Tv. Finally, the Boltzmann equation for the evolution of T v in time due to e + e~ «-> v^v^ can be written as 



dT 2G 2 f°° f e 

r~vT / de 9 el / (fe iJ (e„, e p , e) {(1 - F v )(l - T 9 )T e -T e + - T v Ty(\ - T e -)(A - T e +)} , 
2 *r Jo Jo 



dt (27 



(3.22) 



where e = e v +£„ and H (e„, e„, e) is given in eq. (|B5| ). In solving eq. (|3.22| ), JF P must be evolved simultaneously with 
T v . To do so, in addition to making the appropriate changes to the vector and axial- vector coupling constants, V 
and A, one needs to in tegrate over e„ instead of e p . Note that the electron and positron distribution functions appear 
explicitly in eq. ( 3.22 ). We take these distributions to be Fermi-Dirac at temperature T and with r\ e determined by 
T, p, and Y e . 

Equation ( |3.1§| ) may also be used to find the total volumetric v^v^ pair spectrum by replacing e v in the numerator 
with e. Ignoring neutrino blocking in the final state one can show that pOJ 



2.09 x 10^ 



( T \ 



VMeV 



J f(Ve) ergs cm 3 s \ 



(3.23) 



where 



and 



f(Ve) 



F 4 ( Ve )F 3 (- Ve ) + Fj{-ri e )F 3 {7) e ) 

2^4(0)^(0) 



(3.24) 



F n (y) 



e x-y + l 



■ dx 



(3.25) 



are the Fermi integrals. 



D. Nucleon-Nucleon Bremsstrahlung 



The importance of nucleon-nucleon bremsstrahlung in late-time neutron star cooling has been acknowledged for 



some time [113 14 1. Recently, however, this process has received more attention as a contributor of v u v u pairs and as 



an energy transport mechanism in both core-collapse supernova and nascent neutron star evolution |7|, |12| , pil|]22| | . The 
contribution from nucleon-nucleon bremsstrahlung is a composite of neutron-neutron (mi), proton- proton (pp), and 
neutron-proton (np) bremsstrahlung. Fermi's Golden Rule for the total volumetric emissivity of single v^s due to nn, 
pp, or np bremsstrahlung, including f„ and P„ blocking in the final state, is given by 



Q 



where 



d 3 Pi 



n (27r)3 



d 3 q v d 3 q v 
(27r) 3 2e, (2Tr) 3 2e p 



lu (sEl-^l 2 ) (2vr) 4 <5 4 (P)S[^] 



(3.26) 



(3.27) 



The product of differential phase space factors in eq. ( 3.26| ) includes a term for each of the four nucleon s inv olved 
in the process; 1 and 2 denote initial-state nucleons whereas 3 and 4 denote final-state nucleons. In eq. ( |3.26 ), s is 
a symmetry factor for identical initial-state fermions, q v is the neutrino three-momentum, e„ is the neutrino energy, 
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and the four-momentum conserving delta function is explicit. In a one-pion exchange model for the nucleon-nucleon 
interaction, the spin-summed matrix element can be approximated by p^,|21| 



\M\< 



V ma 



k 2 



1 {e v £„ - q v ■ kq p ■ k) 



(3.28) 



where e = e v + £p, k is the magnitude of the nucleon momentum transfer, gj± — —1.26, / ~ 1 is the pion-nucleon 
coupling, and m w is the mass of the pion. In order to make the 18-dimensional phase-space integration in eq. ( 3.2G| ) 
tractable we assume the quantity in square brackets to be of order unity, but possibly as low as 0.1 |Q. To acknowledge 
our ignorance, we introduce the factor, £, and assume these momentum terms are constant. Furthermore, we neglect 
the momentum of the neutrinos relative to the momentum of the nucleons. We are left with a simple, but general, 
form for the bremsstrahlung matrix element: 



J2\M\ 2 ^AC 



(3.29) 



where A — 64G 2 g\f 4 / 'rn\ . In the case of nn or pp bremsstrahlung, as appropriate for identical particles in the initial 



state, the symmetry factor (s) in eq. ( |3.26 ) is 1/4. Such a symmetry factor does not enter for the mixed-nucleon 
process, np, which is still further enhanc ed b y the fact that a charged pion mediates the nucleon exchange 
This increases the matrix element in eq. (3.29) by a factor of 7/3 in the degenerate nucleon limit and ~ 5/2 in the 
non-degenerate limit [2lf| . Considering the already substantial simplifications made by choosing not to handle the 
momentum terms directly, we will adopt the more conservative 4 x (7/3) enhancement for the np matrix element. The 
total volumetric emission rate combining all processes is just Qtot = Qnn + Qpp + Qnp- What remains is to reduce 
eq. ( |3.26 ) to a useful expression in asymmetric matter and at arbitrary neutron and proton degeneracy. 

Following ref. (|l), we define new momenta, p± = (p\ ±p 2 )/2 and p 3c ^ c = 753,4 — p+, new direction cosines, 71 = 
p+ ■ p- /\p+\\p-\ and 7 C = p + ■pz c /\p+\\pzc\i and let Ui = p 2 /2mT. Furthermore, we note that d 3 p\d 3 p2 = 8 d 3 p + d 3 p-. 
Using the three-momentum conserving delta function, we can do the d 3 p*4 integral trivially. Rewriting eq. ( 3.26| ) with 
these definitions, we find that 



Q = 2AsC(2toT) 9/2 (2tt)- 9 / de v e 3 v \ de p du_ du+ dw 3c d 7i d 7c (e p /e) 2 (u_u+u 3c ) 1/2 6(E) H[JF], 



where 



5(E) =6( J2 * - e) = 6(2T(u_ - u 3c - e/2T)). 



(3.30) 



(3.31) 



i=l 



cosines, the dimensionless momenta (lit), and the initial-state nucleon degeneracy factors 7712 = 



The nucleon distribution functions in the term ^[J-j in eq. (3.30) have been rewritten in terms of the new direction 

(3.32) 



.F1 = 



e 



and 



To = 



where a' 12 = 01,2/2 = \(u^ 



2 cosh(a' 1 + 6' 7 i) "2 cosh(a^ - 6' 7 i) ' 

u_ — 771,2) and b' = 6/2 = (u+u-) 1 / 2 . Furthermore, 

e (c' 1 +d'7 c ) e (c2-d'7c) 



2cosh(ci + d'jc) 



and 



(1 - ^4) 



2 cosh(c2 — d'7 c ) 



(3.33) 



where Cj 2 = c i,2/2 = \(u+ + u 3c — 771,2) and d' = d/2 = (u+usc) 1 / 2 . T v and JF P , in contrast with the nucleon 
distribution functions, are independent of angle; for a given set of thermodynamic conditions, they remain functions 
of energy alone. While non-trivial, the integrations over 71 and 7 C can be performed. For example, the result for the 
71 integration is of the form 



1 



2^B(B + 1) 



: In 



(B - (1 + 2B)£ 2 + 2£Vb(B + 1)(£ 2 - 1) 



(3.34) 



where B = sinh 2 a' and £ = cosh b' 71. With a proper evaluation of the integration limits and some algebra one can 
rewrite this result as 
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2 sinh a' cosh a 1 



In 



( 1 + cosh a cosh b + sinh a sinh b) 



( 1 + cosh a cosh 6 — sinh a sinh fe) 



(3.35) 



Similar operations yield a result for the j c integral in terms of c and d. In addition, eq. (3-31) can be used to eliminate 
the integral over w_. Collectively, these manipulations reveal that the differential bremsstrahlung emissivity at 
arbitrary neutron and proton degeneracy is simply a three-dimensional integral over u + , u 3c , and e p : 



dQ_ 

de v 



Ks( (1 - T v ) e% J de v du+ du 3c (e w /ef u" 1 / V 6 / 2 $(e, u+,u 3c ) (1 - Tp), 



where 



(3.36) 



(3.37) 



$(e,u+,u 3c ) = sinh X (/)ln 
x sinh _1 (j)ln 



1 + cosh(e^ 
I + cosh(e_ 
I + cosh(/i4 



cosh(/) + cosh(g + ) 
cosh(/) -I- cosh(g_) 
cosh(j) + cosh(fc+) 



1 + cosh(ft,_) J \cosh(j) + cosh(fc_) 



(3.38) 



and 



, 1/2 , 1/2-.2 

e± = (u + ' ±u_L ) - t]2 

I = u + + u- - rjx/2 - 772/2 
g ± = ±2(u + u-) 1 ' 2 - 771/2 + 772/2 



j = u + + u 3c - 771/2 - 7/2/2 
fct = ±2(7i+u 3c ) 1/2 - rn/2 + 772/2 



1/2x2 



>]2 



(3.39) 



Though u_ has been integrated out via the energy-conserving delta function, it appears here in an attempt to make 
this expression more compact and should be read as u_ = u 3c + e/2T. Importantly, if 771 = 772 the right-hand term 
within both logarithmic terms in $(e,M+,U3 C ) becomes unity. 

Using eq. (3.21), we can easily obtain the contribution to the Boltzmann equation due to nucleon-nucleon 
bremsstrahlung for arbitrary nucleon degeneracy, in asymmetric matter, and including the full nucleon and neutrino 
Pauli blocking terms. We find that 



jv = K' sC, J dep du + du 3c 



-1/2-06/2 



$(e,u+,u 3c ) (1 - T v ) 



(3.40) 



where K' — [(2ir) 3 /4ir]K . The nucleon phase-space integrations above are identical in form for the Vnpn absorption 
process, v^v^nn — > nn. In this case, then, the primed energies are now associated with nucleons 1 and 2 in the above 
manipulations and the incident nucleons (3 and 4) have unprimed energies. If we take the form derived above for the 
nucleon phase-space terms, the absorption channel (Xv) must then contain a factor of e^ e . In addition, the blocking 
term, (1 — J-„), becomes J-p. The Boltzmann equation for the evolution of T v in time is then, 



\ dT v 
c dt 



= K's( J de„ du+ du 3c ( £p /e) 2 u + 1/2 e~^ 2 $(e,u+,u 3c ) {(1 - F V ){1 - F„) - F V F„ e 0e } 



(3.41) 



For the neutron-neutron (nn) or proton-proton (pp) bremsstrahlung contribution, we simply set s = 1/4 in eq. (3.41) 
and use 7/1 = 772 = r\ n or 771 = 7/2 = ?7p, respectively. For the mixed nucleon (np) bremsstrahlung we set s = 1, multiply 
eq. ( j3~4l| ) by 7/3, and set 771 = r/ n and r/ 2 = r) v . While eqs. fl3.36p and ( |3.41| ) may not appear symmetric in 771 and 
772 the logarithmic terms conspire to ensure that the rates for both np and pn bremsstrahlung are identical, as they 
should be. That is, it makes no difference whether we set 77™ or r\ v equal to 771 or 772. 

Tv must be evolved simultaneously with T v . In this case, however, the 



Just as in § III C , in considering e + e 



situation is simpler. Suppose we wish to compare electron-positron annihilation with nucleon-nucleon bremsstrahlung 
by starting at t — with T v = T v = over all energies. We then solve eq. (3.22) and its T v counterpart at each 



9 



timestep and at each energy. For e + e~ annihilation, T v and T„ will evolve differently; they will be visibly different at 
each timestep, because of the weighting of the vector and axial-vector coupling constants which appear in the matrix 
element. In contrast, eq. (3.41) for bremsstrahlung must be solved only once. Since there is no difference in weighting 
between and v^, we can set Ty = T v at every energy, at every timestep, as long as JF P = T v at t = 0. Of course, if 
we wish to consider T„ ^ T v initially, the two distributions would need to be evolved separately and simultaneously, 
coupled through the blocking and source terms on the right-hand side of the Boltzmann equation. 



Equation ( 3.26| ) can also be used to find the total volumetric v^v^ pair emissivity. To facilitate this we replace e v 



with e and insert J S(e— (e„ + e^))de. Assuming the neutrinos are radiated isotropically, we can use this delta function 
to do the integral over d 3 qp and leave the total rate in terms of an integral over e v from zero to e and another over e 
from zero to infinity. Momentarily ignoring neutrino blocking in the final state, the former can be integrated easily. 
Making the same momentum, angle, and nucleon distribution function substitutions we used in deriving the single 
spectrum we can reduce the pair spectrum to an integral over u + , u% c , and q = e/2T. We find that 



Ds£T 8 ' 5 / dq dusc du+ q A e 



-1/2 



$(e,it+,u 3c 



(3.42) 



where 



D = 



8 G 2 g\ / / 



15 v^TT 9 



,9/2 



(3.43) 



(3.3E 



Note that eq. ( 3.42| ) allows us to easily calculate the pair differential 
For Q"" P(< and Q^ P(J , s = 1/4. As with the single spectrum, for 
multiply eq. (p§) by 7/3 and set a = 1. Finally, Q*f* =QTv u 



and $(e,u+,U3 C ) is defined in eq. 
volumetric emissivity (dQ u u /de) 



1. The N on- Degenerate Nucleon Limit 



In the non-degenerate nucleon limit, the term JFiJ-^l — .7"3)(1 — F4) reduces to e Vl e' l2 e~ 2tyU++u ^ jl2| which is 
independent of angle. This tremendous simplification allows for easy integration over u + and u 3c in eqs. (3.36), 
( |3.41 ), and (3.42). The total volumetric emissivity of a single v^v^ pair in this limit, ignoring and blocking in 
the final state, is |T^] 



1.04 x lO^CpCpu) 2 



5.5 



/ V , | 

ergs cm s . 

MeV/ 8 



(3.44) 



For nn and pp bremsstrahlung, X is the number fraction of neutrons (X n ) or protons {X p ), respectively. For the 
mixed-nucleon process (np), X 2 becomes (28/3)X n X p . Fi g ure |l| compares the non-degenerate nucleon limit (eq. |3.44| ) 
with the arbitrary nucleon degeneracy generalization (eq. [3.42| ) in the case of neutron-neutron (nn) bremsstrahlung, 
as a function of the neutron degeneracy rj n = fi n /l '. The filled square shows the degenerate limit obtained by ref. 
fuj . Note that at r/ n ~ 0, the fractional difference between the two is just ~ 12%. At realistic neutron degeneracies 
within the core (rj n ~ 2), this difference approaches 30%. 

The single differential emissivity can be written in terms of the pair emissivity |l2| ] : 



dQ_ 

de v 



C 



Q"n v » \ 3 



T 4 



-1q v x 



(x 2 -xfl 2 dx 



-Ki(q)(q- q v ) 2 dq 



(3.45) 



where C = 2310/2048 ~ 1.128, q u — £„/2T, q — e/2T, and K\ is the standard modified Bessel function of imaginary 
argument. A useful fit to eq. ( 3.45| ), good to better than 3% over the full range of relevant neutrino energies is [[l2| 



f2 „ 0.234 

de v T 



/ f \2.4 



-l.le„/T 



(3.46) 



Using eq. ( [3.21 ), we obtain the contribution to the Boltzmann equation including Pauli blocking of and v p neutrinos 
in the final state: 
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= CstJ^de* (y) e"^ 2 {(1 - ^)(1 - T„) - T v T„e^} . (3.47) 



dt 

where 



In obtaining eq. (3.4S), we have used the thermodynamic identity in the non-degenerate limit, 

where n is the number density of nucleons considered and i is 1 or 2 for neutrons or protons, depending on which 
nucleon bremsstrahlung process is considered. 



IV. RESULTS 



The numerical algorithm we have developed accepts arbitrary initial and phase-space distributions. Using 
the scattering formalism developed in the previous section, we evolve two initial distribution functions: (1) a broad 
Gaussian in energy centered at 40 MeV with a maximum of T v — 0.80 and a full-width at half-maximum of ~28.6 
MeV, and (2) a Fermi-Dirac distribution at a temperature 2 x the temperature of the surrounding matter and with 
zero chemical potential. While the former is unphysical in the context of supernova calculations, it illustrates the 
effects of blocking on both the average energy transfer and the rates for each scattering process. Furthermore, its 
evolution is more dynamic than the Fermi-Dirac distribution. As a result, the way in which the distribution is spread 
and shifted in time is more apparent. The essential differences between the two processes are then more easily gleaned. 
The latter initial distribution is motivated by consideration of the environment within the core of a supernova. The 
and distribution functions, having been generated as pairs via e + e~ <-> v^v^ and nucleon-nucleon bremsstrahlung 
should have approximately zero chemical potential. Furthermore, even in the dense core, the v^s will diffuse outward 
in radius and, hence, from higher to lower temperatures. By starting with a Fermi-Dirac distribution at twice the 
temperature of the matter at that radius, we learn more about how equilibration might effect the emergent spectrum 
in an actual collapse or protoneutron star cooling calculation. 

For the production and emission processes, we start with zero neutrino occupancy and let each build to an equi- 
librium distribution of v^s and v^s. As a check to the calculation, the asymptotic distribution should be Fermi-Dirac 
at the temperature of the ambient matter with zero neutrino chemical potential. Throughout these simulations, we 
take the factor £ in eq. ( |3.4l| ) for nucleon-nucleon bremsstrahlung to be 0.5. (This factor represents our ignorance of 
the importance of the nucleon momentum transfer terms.) 

We repeat these calculations for four temperature, density, and composition points (StarA, StarB, StarC, and 
StarD) taken from the one-dimensional collapse calculation profile, Star corresponding to four radii below the 
shock (~ 80 km). Roughly, these points have densities 10 14 , 10 13 , 10 12 , and 10 11 g cm -3 . The actual numbers are 
shown in Table |. 



A. Scattering 



Figures ^] and |^ show the evolution of a Gaussian distribution at t — to an equilibrium Fermi-Dirac distribution at 
the temperature of the surrounding matter due to f^-neutron (y^n) and ^-electron (i/ M e~) scattering, respectively. 
The equilibrium distribution has a non-zero neutrino chemical potential set by the initial total number of j/^s, which 
is conserved to better than .001% throughout the calculation. Multiple curves on each plot show snapshots of T v in 
time from t — to 1000 microseconds (us). Both calculations were carried out at the thermodynamic point StarB 
whose characteristics are shown in Table [|. StarB is indicative of the core of a supernova, a region of moderate to high 
temperatures (T ~ 15 MeV) and densities of ~ 10 13 g cm~ 3 . These two figures illustrate the fundamental differences 
between z/ M e _ and v^n scattering as thermalization processes. Curve A in Fig. ^|and curve C in Fig. ^ indicate that 
at high energies (e„ > 30 MeV) v^n scattering is a much more effective thermalization mechanism. At e v ~ 40 
MeV both curves show the distribution is within ~30% of equilibrium. Importantly, however, curve A is at 0.33 /is 
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for v^n scattering whereas curve C is at 3.30 jj,s for v^e" scattering. Curve C, in Fig. || for v^n scattering, also 
at t = 3.30 fis, shows that above ^25 MeV the distribution has almost equilibrated. For v^eT scattering, similar 
evolution at high neutrino energies takes approximately 25 /xs. These simple estimates reveal that v^n scattering is 
about 10 times faster than v^e~ scattering at equilibrating v^s with energies greater than approximately 25 MeV. 

This situation is reversed at low s^s. Comparing curve E at t = 33.0 fis in both Fig. || and Fig. ||, we can see that 
at ~ 10 MeV both distributions have filled to approximately the same percentage of the asymptotic, equilibrium T v . 
However, below e u ~8 MeV, v^n scattering has not filled T v to the extent v^e" scattering has. In fact, the rate at 
which these low energy states are filled by v^n scattering is very low; the energy transfer (ui) is much smaller than 
the incident energy. In this regime, the Fokker-Planck approximation for the time evolution of T v in energy space 
may be applicable. In marked contrast, Fig. || indicates how effective i^e - scattering is at filling the lowest e v states. 
Curves F from Figs. || and pi taken at 1000 [is, show that though the distribution has reached equilibrium via v^e~ 
scattering, for v^n scattering the very lowest energy states remain unfilled. For each of the four points in the Star 
profile we consider, v^n scattering dominates at high energies (> 20 MeV), whereas v^e~ scattering dominates at low 
i/fj, energies (< 10 MeV) and particularly for e v < 3 MeV. 

Figures || and |^ depict the evolution of T v via v^n and v^e~ scattering, respectively, for an initial Fermi-Dirac 
distribution at 2x the temperature of the surrounding neutrons and electrons and with zero neutrino chemical po- 
tential. This calculation was carried out at StarC (see Table |), which is representative of the outer core, in the 
semi-transparent regime, where the neutrinos begin to decouple from the matter (near the neutrinosphere). The same 
systematics highlighted in the discussion of the evolution of the initial Gaussian distribution for StarB are borne out 
in these figures. Curves A and B on both plots, denoting 0.10 and 0.33 milliseconds (ms) of elapsed time, respectively, 
confirm that above e„ ~ 15 MeV v^n scattering dominates thermalization. 

Figures [| and ^ show (u>)i n and (u>)outi as defined in eqs. ( |2.6| ) and (2.7), for v^n scattering and v^e~ scattering, 
respectively. The separate curves portray the evolution in time of the thermal average energy transfers as the 
distributions evolve to equilibrium (cf. Figs. [| and ||) . As one would expect from kinematic arguments, the magnitudes 
of both (u>)i n and {u>) ou t for v^n scattering are much less than those for v^e~ scattering. Though the energy transfers 
are much smaller, even at the highest energies, v^n scattering still dominates v^e~ scattering in thermalizing the 
distribution because the rate for scattering is so much larger. At low neutrino energies, however, both average energy 
transfers for neutron scattering go to zero, whereas they approach large negative values (~ —20 MeV) for electron 
scattering. At these low energies, the fact that the rate for v^n scattering is larger than for ^ M e~ scattering fails to 
compensate for the vanishing energy transfer. For example, at e v — 3 MeV and t = 33 ms, the energy transfer for 
v^e" scattering is more than 100 times that for v^n scattering. 

In order to fold in information about both the rate of scattering and the average thermal energy transfer, we plot 
Tr> and Te (eqs. 2.8 and 2.S) in Fig. || for all four points considered in the Star profile. We show here a snapshot of 



the rates for both scattering processes for a Fermi-Dirac distribution initially at twice the local matter temperature, 
with zero neutrino chemical potential. Note that the spikes in indicate the neutrino energy at which {ui) ou t — 
(cf. Figs. | and |). In general, we find that as e v — > 0, and Te go to zero for ^ M -neutron scattering, whereas 
Td approaches a constant and Te gets very large for v^e~ scattering |l5| . This is a consequence of the fact that, 
regardless of T v , Uo) ou t — * for v^n scattering as e v — * 0, as shown in Fig. ^|. For v^e~ scattering the situation is 
different. As Fig. |7] reveals, {uj) out approaches ~ —20 MeV at e u — 0. As expected from our analysis of the evolution 
of T v in Figs. [|and ||j a t approximately 40 MeV the thermalization rate for v^n scattering for StarB is about an order 
of magnitude greater than that for v^e~ scattering. Specifically, the Trj's cross at ~15 MeV, whereas the Te's cross 
at ^20 MeV. Below these energies, both v^n rates drop off precipitously as a consequence of the fact that {u))out — ► 0. 
Below e v ~ 5 MeV, the thermalization rate for v^e~ scattering dominates by 2-5 orders of magnitude. As evidenced 
by the other panels in Fig. ||, this same trend holds in the other regions of the stellar profile. In general, the rates drop 
over the whole energy range for both processes as the density and temperature decrease, but the same systematics 
hold. In fact, for StarA, StarC, and StarD the Te and crossing points for both processes are lower than those 
for StarB. As a result of the higher temperature at this radius (T ~ 14.5 MeV) v^e" scattering is important in 
thermalizing slightly higher energy neutrinos than at the other radii. For StarC and StarD, specifically, both rates 
cross at neutrino energies less than 12 MeV. 

These results demonstrate that ^-nucleon scattering is an important thermalization process from the dense core 
through the semi-transparent regime for v^s with energies greater than approximately 15 MeV. The addition of this 
energy transfer mechanism implies that the u^s stay energetically coupled to the surrounding matter longer than has 
been previously estimated |23[ |. We can approximate the radius at which the f M s energetically decouple from the 
matter (the Ensphere) p3| by observing when the diffusion timescale is approximately equal to the equilibration 
timescale given by T^ 1 = td, as defined in eq. (2.8). Using this crude approximation we find that by including 
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i^-nucleon energy transfer the E^-sphere is pushed outward in radius by approximately 3 kilometers. This difference 
in radius corresponds to a 1-2 MeV drop in the matter temperature in the model Star. The average energy of the 
emergent spectrum is roughly correlated with the local matter temperature of the ^-sphere. Therefore, we conclude 
that i/^-nucleon energy transfer in full transport calculations will likely soften the emergent v n spectrum. 



B. Emission and Absorption 

Figured shows the total integrated volumetric emissivity as a function of radius in the model Star for nucleon-nucleon 



brcmsstrahlung in the non-degenerate nucleon limit (eq. 3.44), its generalization for arbitrary nucleon degeneracy (eq. 
3.42), and the emissivity for e + e~ annihilation (eq. 3.23| ). Note that not one of these expressions contains neutrino 



blocking terms and that the general bremsstrahlung rate crosses that for e + e _ annihilation at ~ 23 kilometers where 
p ~ 6x 10 12 g cm -3 , T ~ 11 MeV, and Y e ~ 0.13. While this plot gives a general idea of where e + e~ annihilation should 
begin to compete with nucleon-nucleon bremsstrahlung, it fails to include the differential nature of the production in 
energy. In addition, it does not include absorption or blocking effects, which quantitatively alter the relative strength 
of the emission. 

To begin to understand the import of these terms and the character of each pair production process, we include 
Figs. |l0| and 11, which show the time evolution of T v via nucleon-nucleon bremsstrahlung and electron-positron 



annihilation, respectively, for the point StarC, initialized with zero v n and Sy, phase-space occupancies. The final 
equilibrium distribution is Fcrmi-Dirac at the temperature of the surrounding matter, with zero neutrino chemical 
potential. Comparing curve C on both graphs, which marks 10.0 milliseconds (ms) of elapsed time, one can see that 
bremsstrahlung dominates production below ~ 15 MeV. Indeed, bremsstrahlung overshoots its equilibrium distribu- 
tion at energies below 10 MeV before finally filling the higher e v states. In contrast, electron-positron annihilation 
fills the higher states first and moves slowly toward the low-lying neutrino energies, taking a factor of 10 more time 
at this thermodynamic point to reach equilibrium. 



In Figs. [12| and |13|, we plot Ti„ and T out , as defined in eqs. (2.4) and (2.5), for both production processes at 
the point StarB. As one would predict from our simple observations of the time evolution of T v , the bremsstrahlung 
rates are much faster (~2 orders of magnitude) than the e + e~ annihilation rates at low neutrino energies. At StarC, 
e + e~ annihilation competes with bremsstrahlung above e v ~ 15 MeV. For StarB, however, at a matter density an 
order of magnitude greater than that for StarC, the energy at which nucleon-nucleon bremsstrahlung becomes more 
important than e + e~ annihilation is ~ 60 MeV. In this regime, where T ~ 12—14 MeV and p ~ 10 13 g cm" 3 , 
we find that bremsstrahlung dominates neutrino pair-production via electron-positron annihilation. A close look at 
the evolution of the total thermal average neutrino energy {{e v )) reveals that T v reaches its asymptotic equilibrium 
distribution via nucleon-nucleon bremsstrahlung in ~ 1 ms. Electron-positron annihilation takes ~ 50 ms to fill all 
but the very lowest energy states. This trend continues as the matter becomes more dense. For StarA, well beneath 
the neutrinospheres at p ~ 10 14 g cm -3 , the rates for bremsstrahlung and electron-positron annihilation never cross. 
In fact, the former produces an equilibrium Fermi sea of v^'s in ~ 50 /is, whereas the latter takes ~ 10 3 seconds. 
This difference of 8 orders of magnitude in timescale, however, is a bit misleading. Similar to v^n scattering, e + e~ 
annihilation has trouble filling only the very lowest neutrino energy states. In actuality, at the highest energies, both 
Ti n and T out for e + e~ annihilation come within 3-4 orders of magnitude of the rates for bremsstrahlung at the same 
energy. Still, the difference is striking. As the temperature drops from StarB (14 MeV) to StarA (10 MeV) and the 
density increases by an order of magnitude, rj e goes from 3.79 to 15.75. Consequently, Pauli blocking of electrons in 
the final state suppresses the process v u v u — > e + e~ , and the phase-space density of positrons is depleted to such an 
extent that e + e _ — > v u v u is suppressed as well. We conclude that beneath the neutrinospheres and specifically for 
p 10 13 g cm -3 , nucleon-nucleon bremsstrahlung is the primary and dominant v u v^ source. Near the neutrinosphere, 
within the gain region and behind the shock, between 30 km and 60 km at p ~ 10 12 g cm -3 and T ~ 6 — 8 MeV, 
bremsstrahlung competes with e + e~ annihilation at all neutrino energies and is the primary production process for 
the low-lying e v and e p states. 

The addition of nucleon-nucleon bremsstrahlung will have quantitative implications for the and v T emergent 
spectra. Specifically, they should be softer and brighter. Burrows et al. |l2j confirm this with their study of static 
supernova and protoneutron star atmospheres, having included nucleon-nucleon bremsstrahlung in the non-degenerate 
limit. In addition to observing a systematic softening, they also find that the spectrum is a factor of 2 more luminous 
at e v = 10 MeV. 
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V. SUMMARY AND CONCLUSIONS 



Our results for equilibration via f^-electron scattering and f^-nucleon scattering indicate that the latter competes 
with or dominates the former as a thermalizer for neutrino energies > 10 MeV for p > 1 x 10 11 g cm" 3 at all 
temperatures. At neutrino energies > 30 MeV the difference at all densities and temperatures is approximately an 
order of magnitude. For the production and absorption processes, we find that nucleon-nucleon bremsstrahlung, 
at the average energy of an equilibrium Fermi-Dirac distribution at the local temperature, is 5 and 2 orders of 
magnitude faster than e+e" annihilation at StarA (T ~ 10 MeV, p ~ 10 14 g cm" 3 ) and StarB (T ~ 15 MeV, 
p ~ 10 13 g cm -3 ), respectively. Only for p ~ 10 12 g cm" 3 and T ~ 6 MeV does e + e~ «-> v^v^ begin to compete 
with bremsstrahlung at all energies. We conclude from this study that the emergent and v T spectrum is (1) 
brighter and (2) softer than previously estimated. The former results from the inclusion of the new pair emission 
process, nucleon-nucleon bremsstrahlung. The latter is a consequence of both the increased energy coupling between 
the nuclear and neutrino fluids through ^-nucleoli scattering and the fact that bremsstrahlung dominates e + e~ 
annihilation near the neutrinospheres at the lowest neutrino energies. While the full transport problem, including 
i/^-nucleon scattering energy redistribution and nucleon-nucleon bremsstrahlung, must be solved in order to delineate 
precisely what consequences these processes have for the emergent spectrum, these calculations demonstrate that 
they should not be omitted. 
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APPENDIX A: NEUTRINO-ELECTRON SCATTERING 



Each of the retarded polarization functions in eqs. (3.15-3.17) can be written in terms of one-dimensional integrals 
over electron energy (e e ), which we label I n j|; 
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(A2) 



hnn*t{q,u) = _/ 0) (A3) 

and 

ImlT^ (q, u) = Ko + 2I X ]. ( A4) 

The authors of [0 were able to express the 7„'s in terms of polylogarithmic integrals such that 

la =Tz(l-^\, (A5) 



h = T 2 z 
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zT 



(A6) 



and 
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/ t-< I 2 _l 71-2 , z2 _lo^ 3 e-6 , ejg x 



(A7) 



where 7y e = /i e /T is the electron degeneracy, z = /3w, is the energy transfer, and 



In eqs. (A5-A7), the £ n 's are differences between polylogarithmic integrals; = Li n (— ai) — Li(— 012), where 



Jo £ 



(A8) 



(A9) 



and Lii(a;) = ln(l — a;). The arguments necessary for computing the integrals are a.\ = cxp[/3(e_ + u>) — rj e ] and 
a 2 = exp(/?e_ - r/ e ). 



APPENDIX B: ELECTRON-POSITRON ANNIHILATION 



The production kernel is defined by 



R p (s v , Ev, cost 



d 3 p d 3 p' 



2£„£ p 7 (27r) 3 2e (27 

The differential production spectrum for final state v^s can then be written as |2| 

^■ = ( 1 -^)(^6 £2d£ -/ ^'^ d^B?{E Vl e v ,cos6){l-Tv) 



(Bl) 



(B2) 



where df2 is the differential solid angle for the neutrino, \i! = cos 6' is the cosine of the angular coordinate, and 
<f> is the azimuthal angle between and v^. Expanding the production kernel in a Legendre series in the scattering 
angle, cos6> = ^ + [(1 - /i 2 )(l - // 2 )] 1/2 cos (f>, 



R p (e„,e„cos0) = \^ 21 + Vtf&i ^W™ 3 e ) ~ e p ) + \${{e v ,e v ) cos 6, 



(B3) 



$0, in eqs. (3.20) and (BS), is given by 



q2 



(B4) 



(B5) 



where !F e + is a function of e'{= e v + £y — e) and 

H (e v ,£v,e) = (V + A) 2 Jq(e u ,e p , e) + (V — A) 2 J f7 (e I/ , e p , e) . 
Each Jo in eq. (B5) is a polynomial in e„, Sp, and e of dimension [energy]. They are related to each other by j^] 

Jo(e„, En, e) = J r/ (£p, s v , e) . (B6) 



Both Jq and Jq 1 can be found in ref. ^ (correcting for the typo in their eq. C67). From eqs. (3.20) and (B5) we see 
that the differences between the spectra for v^s and v^s for a given temperature and electron degeneracy (r} e ) arise 
solely from the relative weighting constants (V + A) 2 and (V — A) 2 in eq. (B5) for Jq and Jq 7 , respectively. Indeed, 
in this approximation the same can be said for the difference in the spectrum between u e and neutrinos. 
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TABLE I. Radius (R), temperature (T), density (p), lepton fraction (Y e ), and degeneracy factors (ij = (J./T) for neutrons, 
protons, and electrons for four points from the model, Star, a one- dimensional core-collapse calculation evolved through collapse 



for 0.24 seconds 


. At this point 


in the core evolution, 


the shock is 


at about 80 kilomet 


ers Burrows, 


Hayes, and 


Fryxell (1995). 


Label 


R (km) 


P (g cm" 3 ) 


T(MeV) 


Ye 




Vp 


Ve 


StarA 


10.75 


1.281 x 10 14 


10.56 


0.2752 


2.37 


0.70 


15.75 


StarB 


18.75 


1.023 x 10 13 


14.51 


0.2021 


-1.62 


-3.04 


3.79 


StarC 


34.75 


1.082 x 10 12 


6.139 


0.0907 


-2.48 


-4.81 


3.03 


StarD 


49.75 


1.071 x 10 11 


4.527 


0.1671 


-4.45 


-6.06 


1.93 
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FIG. 1. The total volumetric emissivity 
non-degenerate neutron limit (solid line, eq. 
MeV, Y e = 0.0 



due to neutron-neutron bremsstrahlung (Q„ n ) in 



ergs 
eq. 



3.44) and at arbitrary nucleon degeneracy (dashed line 
and for a range of densities from 5 x 10 10 g cm -3 to nuclear density (~ 2.68 x 10 14 g cm 



3.42) for 



the 
= 6 



3 ). The filled box 



denotes the degenerate neutron limit obtained by Flowers, Sutherland, and Bond (1975). 
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FIG. 2. The time evolution via ^-neutron scattering of the neutrino distribution function (J- v ) for an initial Gaussian 
distribution centered on 40 MeV, for the thermodynamic characteristics specified by StarB in Table @. The curves show the 
distribution at snapshots in time: (A) t — 0.33 /is, (B) t = 1.00 [is, (C) t — 3.30 /is, (D) t = 10.0 (is, (E) t = 33.0 ^s, and 
(F) t = 1000 ^s. The solid dots denote an equilibrium Fermi-Dirac distribution at the temperature of the surrounding thermal 
bath with a neutrino chemical potential [i v ~ 2.32T set by the initial neutrino number density. 
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FIG. 4. The time evolution via i/ M -neutron scattering of the neutrino distribution function {T u ) for an initial Fermi-Dirac 
distribution at 2x the ambient temperature, for the thermodynamic characteristics specified by StarC in Table [[]. The curves 
show the distribution at snapshots in time: (A) t = 0.10 milliseconds (ms), (B) t = 0.33 ms, (C) t — 1.0 ms, (D) t = 3.3 ms, 
and (E) t = 33.0 ms. The solid dots denote an equilibrium Fermi-Dirac distribution at the temperature of the surrounding 
thermal bath with a neutrino chemical potential ~ 2.55T set by the initial neutrino number density. Comparison of this 
plot with Fig. ^ shows that v^n scattering dominates thermalization above e v ~ 10 MeV. 
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FIG. 5. The same as Fig. ^, but for ^-electron scattering. Comparison of this plot with Fig. |4| shows that ^-electron 
scattering dominates thermalization below e„ ~ 10 MeV. 
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FIG. 6. The thermal average energy transfers, (ui)i n and {ui) ou t, defined in eqs. (2.6) and (2.7), respectively, as a function of 
neutrino energy (e v ) for f M -neutron scattering at the thermodynamic point StarC. The curves show snapshots of the average 
energy transfers in time as T v evolves (see Fig. ^). For (a>)j n , (A) t = 0.0033 milliseconds (ms), (B) t = 0.33 ms, and (C) 
t — 33.0 ms. We show (u))out at t = (thin line) and t — 33.0 ms (thick line). Note that in equilibrium (t ~ 33.0 ms) 
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FIG. 8. Td and Te as defined in eqs. (2.8) and (2.S ), respectively, for both v^n and i^e - scattering for an initial Fermi-Dirac 
distribution at StarA, StarB, StarC, and StarD (see Table |j) at a snapshot in time. The spikes in the Td curves are a consequence 
of the fact that {(J) ut — > at those neutrino energies (compare the plot above for StarA with Figs. ^ and [?]). The solid and 
the short-dashed lines in all four plots are Te and Td, respectively, for v^n scattering. The long-dashed and long-short-dashed 
lines are Te and Td, respectively, for f^e - scattering. At all four points in the stellar profile Star, v^n scattering dominates 
v^e~ scattering at energies above 10—20 MeV by approximately an order of magnitude. For StarA, the points where the rates 
for ^-electron and i/ M -neutron scattering cross are at ~ 8 and ~ 13 MeV. For StarB, they lie higher, at ~ 15 and ~ 20 MeV, 
due predominantly to the higher temperature at that radius. For both StarC and StarD, the rates cross at between 7 and 12 
MeV. 
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FIG. 9. The integrated total volumetric emi ssivit y in ergs cm" 3 s _1 for e + e~ annihilation and nucleon-nucleon 
bremsstrahlung in the non-degenerate limit (Eq. 3.44) and at arbitrary degeneracy (Eq. 
the stellar collapse profile, Star. Note that the total bremsstrahlung and e + e~ 



3.36) as a function of radius in 
annihilation emissivities cross at R ~ 23 



kilometers where p ~ 6 x 10 g cm and T ~ 11 MeV. Above this radius, 



dominates. 
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FIG. 10. The time evolution of T v due to nucleon-nucleon bremsstrahlung via eq. (3.41) for the point StarC described in 
Table Q starting with T v — Tv — at all energies. Curves show snapshots of the evolution of T v with time: (A) t = 0.66 
milliseconds (ms), (B) t = 3.30 ms, (C) t = 10.0 ms, (D) t = 33.0 ms, and (F) t = 330 ms. The solid dots denote 
an equilibrium Fermi-Dirac distribution at the temperature of the surrounding thermal bath with zero neutrino chemical 
potential. 
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FIG. 11. The same as Fig. |ly, but for e + e annihilation via eq. (3.22). Curves show snapshots of the evolution of T v with 
time: (C) t = 10.0 milliseconds (ms), (D) t = 33.0 ms, (E) t = 100.0 ms, (F) t = 330.0 ms, (G) t = 1000 ms, and (H) 
t = 3300 ms. 
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FIG. 12. Yin as defined in eq. (2.4) for nucleon-nucleon bremsstrahlung (solid lines) and e + e~ annihilation (dashed lines), 
for the point StarB. Each curve shows a snapshot of Ti„ as T v builds from zero phase-space occupancy at t — 0: for (A) and 
(A') t = 10.0 ps. (B) and (B') denote t = 33.0 /is. (C) marks the equilibrium rate for bremsstrahlung at t — 3.3 milliseconds. 
Curves (D) and (E) mark 100.0 us and 66.0 milliseconds, respectively, for e + e~ annihilation. The latter, marks the e + e~ 
equilibrium rate. Note that the equilibrium rates cross at e v ~ 65 MeV. 
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FIG. 13. The same as Fig. for the same times, but for Tout as defined in eq. (2.5) 
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